Scale free networks from a Hamiltonian dynamics 
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Contrary to many recent models of growing networks, we present a model with fixed number of 
nodes and links, where it is introduced a dynamics favoring the formation of links between nodes 
with degree of connectivity as different as possible. By applying a local rewiring move, the network 
reaches equilibrium states assuming broad degree distributions, which have a power law form in 
an intermediate range of the parameters used. Interestingly, in the same range we find non-trivial 
hierarchical clustering. 



PACS numbers: 05.10.-a,89.75.Hc,05.65+b 

In their theory of random graphs (RG), Erdos and 
Renyi showed that these graphs, composed of N vertices 
(or nodes), connected probabilistically by a set of edges 
have several interesting properties 0. Among them the 
most striking one is the slow rate of growth (~ log N) 
of the diameter of the giant component. This "small 
world" property is very important in connected networks 
represented by single component graphs, since it reflects 
the efficiency of the network for transport or communica- 
tions 0- Over last few years it is becoming increasingly 
evident that most real-world networks have indeed small 
world properties 0, 0], e.g., electronic communication 
networks like Internet 0, World-Wide Web @, social 
networks of acquaintances Q and of collaborations 

On the other hand, some important properties distin- 
guish real-world networks from RG, motivating the rapid 
growth of interest in this field. Many real-world net- 
works have broad nodal degree distributions, P(k) (the 
degree of a node is the number k of links meeting at that 
node) often characterized by a power tail, P(k) ~ fc~ 7 , 
that indicates a scale free (SF) character of the network 
(SFN) pj. Moreover, in real- world networks one observes 
a high degree of clustering, which measures the local cor- 
relations among the links of the network and implies that 
neighbors of a node are more likely to be neighbors |2( 
(this feature has also been associated to the term small 
world 0). The clustering often scales with the degree 
of the relative node. This is connected to a hierarchical 
organization of the network 0, where clustered blocks 
connect to form larger units, and etc. 

To introduce the correlations between nodes that dis- 
tinguish SFN's from RG's, in the last years SFN's have 
been extensively modeled by growing networks in which 
a preferential attachment (PA) rule shapes the nodes 
degree 0, 0, (i.e., each new node is linked to an old 
one with a probability proportional to the degree of the 
old node llCI) . However, biological networks, including 
food webs metabolic networks [r3.ll3| and protein- 
protein interaction networks [l4l ITEj display the features 
listed above, although both the PA and the growing 
process are debatable in these cases. For example, for 
protein-protein networks they have proposed both grow- 
ing network models without PA [16J and models with a 



dominant stationary, asymmetric PA 01 ■ I n food webs, 
where links represent the prey-predator relations, the PA 
and growth are particularly unsuitable to describe the sit- 
uation. Thus, in order to achieve a better understanding 
of the principles shaping a part of the real networks, it is 
worth to spend some efforts to discover dynamics lead- 
ing to non-growing SFN's. It has been already shown 
that models with fixed number of nodes do not require 
linear PA, but SF distributions arise from an algebraic 
PA rule including an exponent which can vary in a wide 
range [l8| . 

In this paper we address the particular problem of 
whether SFN's can arise from mechanisms excluding 
both PA and growth. Recent works 0, I20I l2l| proposed 
theories of networks at equilibrium. In some of these 
cases [gjj, El a SFN can be generated simply by choos- 
ing the desired degree distribution to be SF. On the other 
hand, different works obtained SFN without plugging in 
an a priori degree distribution H |H H E| ■ In the 
spirit of the latter strategy, we propose an example of 
equilibrium network with Hamiltonian that can yield hi- 
erarchical SFN's. The energy function depends on the 
degrees of the nodes and of their neighbors . Since it 
favors connections between nodes with degrees as differ- 
ent as possible, it leads to networks with disassortative 
mixing 26]. Furthermore, the simulation is implemented 
by using a local rewiring rule. Dynamics of this kind 
appear natural for biological networks, which indeed are 
disassortative |l3U26l] . In particular, in food webs we ex- 
pect each species to find not convenient to interact with 
similar (and competing) ones. We notice that in food 
webs both exponential and SFN's are found II], as we 
have in our model. 

We consider a connected network, represented by a sin- 
gle component undirected graph, composed of N nodes 
connected by L undirected links (edges). The network 
topology is uniquely determined by its adjacency ma- 
trix c, such that cy = 1 if nodes i and j are linked, 
and otherwise. The degree of a node i is indicated as 
ki = J2j^i c ij ■ We define an energy associated with a link 
between the i-th and the j-th nodes having the following 




FIG. 1: (a)-(d): Example with N = L = 64, k = 2 at a = 2. Snapshots after £ iterations of the LRM (f), with (a) t = 0, (b) 
i = 4JV, (c) t — 12N and (d) t = N 2 . Blue (darker in gray scale) links meet at the node with the highest degree, (e): Different 
layout of (d). (f): LRM described in the text. 



form: 



min{ki, kj} 
maxjfci, kj} 



(1) 



This equation implies that the energy of a link decreases 
with the difference in nodal degrees at the two ends of the 
link and it contributes no energy when the link connects 
two nodes of same degrees. The Hamiltonian of a network 
configuration Q is then 



(2) 



i<3 



To generate the initial connected network we first add 
Lq = N links to form a graph with the topology of a 
ring (see Fig. da)). The remaining L — Lq links are 
added sequentially and randomly, to connect unlinked 
nodes. The network is evolved by using a local rewiring 
move (LRM), depicted in Fig. 2Jf). A set of three nodes 
is randomly selected. First, a node i is selected with 
probability l/N, secondly the node j which is a neighbor 
of i is selected with probability and finally the node 
I 7^ i which is a neighbor of j is selected with probability 
l/(kj — 1). If cu = 0, a LRM attempts to delete the i-j 
link of the graph Q and introduce the i-l link to obtain 
the new graph Q' with a probability 



P 



-a[H{G') 



(3) 



where a is a tunable parameter. For a — the LRM is 
always accepted [l^. In this case the difference between 
the typical graphs and RG's is reflected in the degree 
distribution. Using a master equation approach |28j | , one 
can show that the degree distribution P{k) indeed decays 
exponentially, as shown numerically in Fig. [21 In the op- 
posite limit of a — > oo, LRM strongly favors connecting 
nodes with degrees as different as possible. As a result, 
usually graphs have several high degree nodes (hubs) con- 
nected to many other mono-degree nodes (leaves) . Notice 
that the LRM cannot split the network in disjoint com- 
ponents. 

The probability |JSJ, the introduction of the Hamilto- 
nian (J2J and the LRM have been chosen for their sim- 
plicity and for their analogy with usual rules of equi- 
librium statistical mechanics, but they do not yield an 



equilibrium distribution of © at an inverse tempera- 
ture a. However, one can see that the above implemen- 
tation gives the canonical distribution of configurations 
with weights Y\^ =1 {ki — 1)! e^ aIi ^\ hence allowing a de- 
scription of the system in terms of ensemble of connected 
networks at the equilibrium given by the Hamiltonian 
H eq (g) = -a -1 J^iLi m ( fc * _ !) ! + H{Q), at tempera- 
ture a^ 1 . Thus, H(Q) can be thought as an interaction 
term from which we expect the arising of complex cor- 
relations in the network. Due to H(Q), the fraction of 
second neighbors of the node i represented by the first 
neighbors of node j and of node I (in the LRM, Fig. [Iff)) 
belongs to the subset of nodes contributing to the energy 
balance in a LRM. The non-trivial build up of the corre- 
lations necessary to obtain SFN's should require 7 < 3, 
meaning an average number of second neighbors ~ N. 
On the contrary, 7 > 3 would forbid the LRM to "feel" 
the global structure of the network, hardly giving a fine 
self-tuning of the network correlations. As shown below, 
indeed we find SFN's with 2 < 7 < 3. 

We apply the described dynamical rules to networks 
composed by N up to 8192 nodes and L = Nk/2, where 
the average degree k is fixed. Thus, we consider sparse 
networks in the limit of large N, where the number 
of links L is much smaller than the maximum number 
L max = N ( N _ x y 2 of possible links in a TV-clique. This 
is motivated by the case of most of the real networks, 
where typically a link between two nodes is an expensive 
or rare object. In our case, we have mainly used k = 4, 
such that L = 2N. 

After a large number r « 2A 2 of LRM attempts on 
the initial configurations with k = 4, we observe a signif- 
icant reorganization of the entire network structure. For 
t <; r iterations of the LRM, equilibrium is reached, as 
indicated by the stable shape of the distribution P Qj Ar(fc) 
of the degrees (see Fig- EJ- The most interesting region 
is 0.8 < a < 4.0, where P a ,N(k) appear as power laws. 
By increasing a, their slopes decrease and their shape 
changes at large fe's, where a shoulder grows, indicat- 
ing the enhanced tendency in the network to form hubs, 
as expected. For a > 4.0 the fraction of hubs is even 
more consistent and the shoulder at high degrees in the 
Pa,N{k) is substituted by a bump (see for instance the 
curve for a — 8 in Fig.|2Jl. Contrary to other models [29| . 
an analysis of the mean degree of the largest hub indi- 
cates that here it does not attract a finite fraction of 
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FIG. 2: Log- log plot of the degree distributions for N = 4096, 
k = 4 and at various q's. 

links, for N — > oo. 

We now focus on the range 0.8 < a < 4.0, where the 
degree distributions are broad and power-law like. To 
support this SF picture we plot Figure in which a 
remarkable feature of this model is evident, namely, the 
self-similar shape of the P a ,N(k) for fixed a. This is 
consistent with a finite size form 

P a . N (k) ~ k~^f (k/N u ^) , (4) 

where f(x) is a scaling function giving a cutoff for suffi- 
ciently large values of x, and f(x) ~ consf for a; — ► 0. In 
order to quantify the SF nature of the networks and to 
support the proposed scaling we perform a finite size 
analysis, first extrapolating the value of the distributions 
slope at fixed a and for N — > oo. In Table [I] we show 
the results. For each a, the value of 7 so obtained is the 
starting point for attempting a data collapse, by using 
a rescaling of the form N lv P(k) vs k/N u . The values 
of 7 that give better collapses (see Table [fl turn out to 
be close to the extrapolated values, such that the whole 
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FIG. 3: Degree distributions for a = 1. Inset: distributions 
rescaled (3) with the 7 and u values quoted in Table 
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FIG. 4: Mean diameter T> and mean clustering coefficient C 
as a function of iV (linear- log and log-log plot, respectively), 
for the same set of a values. 

picture is consistent. In addition, we note that the data 
collapses get worse close to the boundaries of the SF re- 
gion 0.8 < a < 4.0, while outside this range we can not 
make good rescalings. This supports our first, subjective 
delimitation of the SF region (we stress that the k = 4 
case is treated). 

It is interesting to examine the other stylized features 
of networks. In Figure0]we plot the mean diameter T>(N) 
(the maximum of the shortest paths between any nodes 
of a network) as a function of log N for the some represen- 
tative a values. The curves are consistent with a sca ling 
T>(N) ~ A + B\ogN (it seems to be sub-logarithmic [3(j 
for high a). Thus, not surprisingly the small world pic- 
ture is recovered also in our model. 

Given a node i connected with fcj neighbors, if m, is 
the number of links between these neighbors, one can 
quantify the (local) degree of clustering by the clustering 
coefficient C, = mi/L™ ax . The mean clustering coeffi- 
cient C is then the average of the Cj over all nodes of all 
graph realizations at a given a. In Figure 0] we also plot 
C as a function of N, in a log-log scale. The plots are 
compatible with C - N^ ^ for N -> 00, with a rang- 
ing from « 1 to w 0.25 for < a < 6. For the studied 
N's, we notice that the highest clustering is found close 
to a = 4.0, while it is smaller for high and low a. 

Our model also shows the power law dependence 

C{k) - k~ li (5) 

indicating that nodes with few links are typically well 
clustered while hubs hardly are related to high clustering. 
The concept of modularity was introduced to account for 
the hierarchical clustering found in many networks 0, 



TABLE I: Numerical evaluation of the exponents of Eq. @. 



a 


0.8 


1 


2 


3 


4 


7° 


2.9(2) 


2.8(2) 


2.4(1) 


2.3(2) 


2.1(3) 


7 6 


3.0 


2.8 


2.3 


2.2 


2.2 


u b 


0.4 


0.45 


0.55 


0.58 


0.58 



"Values extrapolated for TV — ► 00. 
''Values that give the best data collapse. 
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FIG. 5: Clustering coefficient as a function of the degree, for 
four values of a and N = 8192. 

Isij . In this contexts, networks are built of rather well 
identifiable clusters, which are themselves composed by 
clustered subunits, and so on. While in Figure we see 
that j3 = for a — 0, as for RG, for increasing a a scaling 



(0 takes place in a non negligible interval of k, with a 
(3 raising with a up to values close to j3 = 1. Hence, 
in a present issue on whether = 1 is universal Q, our 
result is in favor of the non-universal character of this 
exponent. 

In summary, we have described a static network model 
with a dynamics favoring networks with disassortative 
mixing and high clustering, where at least three phases 
can be identified by increasing the parameter a, respec- 
tively: exponential, scale free and hub- leaves regimes. 
The scale free regime has a signature of modularity (the 
clustering coefficient scales as a power law of the degree 
with a non-trivial exponent), and the exponent 7 is com- 
prised in the interesting range 2 < 7 < 3. Thus, many 
of the characteristics displayed by real networks, usually 
associated to growing networks with preferential attach- 
ment, can be obtained as well in networks with a fixed 
number of nodes, by using a random rewiring that does 
not require preferential attachment. 
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